IASSNS-HEP-95/62 
Statistics Of The Burst Model At Super-critical Phase 



T. C. Chan**, H. F. CW § , and K. S. Cheng* 
'Department of Physics, University of Hong Kong, Pokfulam Road, Hong Kong 
^School of Natural Sciences, Institute for Advanced Study, Olden Lane, 
Princeton, NJ 08540, USA 
(February 1, 2008) 



Abstract 

We investigate the statistics of a model of type-I X-ray burst [Phys. Rev. E, 
51, 3045 (1995)] in its super-critical phase. The time evolution of the burnable 
clusters, places where fire can pass through, is studied using simple statistical 
arguments. We offer a simple picture for the time evolution of the percentage 
of space covered by burnable clusters. A relation between the time-average 
and the peak percentage of space covered by burnable clusters is also derived. 
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I. INTRODUCTION 



In recent years, spatially extended systems of many degrees of freedom have attracted 
considerable attention. A number of cellular automaton models, for example, Eden model 
0, percolation @, contact process 0, and various self-organized critical models HH have 
been investigated. The forest fire model (FFM) || and some sandpile models with special 
geometries or toppling rules || tell us that the existence of local particle conservation law is 
not a necessary condition for self-organized critical (SOC). Thus, some of the SOC models 
may be used to comprehend the underlying principles of dissipative systems. 

As non-conserving systems are divers in many branches of study, we have proposed a 
simple cellular automaton model for diffusive and dissipative systems @, whose rules are 
similar to that of the FFM. Our model is applicable to a number of phenomena including 
the type-I X-ray bursts ||, and the CO2 gas outburst in some crater lakes ||. Because 
of this, we name our model the "X-ray burst model", or simply the "burst model". We 
studied the burst model by both numerical and semi-analytical means and we found that 
the model exhibits various phases when choosing different set of parameters In spite 

of the similarity with the FFM, the burst model is not self-organized critical. 

In this paper, we focus on the statistical behavior of the burst model at its super- 
critical phase. First, we briefly review the burst model in Section II. Then we discuss the 
relationship of the burst model at its super-critical phase with percolation in Section III. A 
relation between the time-average coverage fraction and the average coverage fraction just 
before burning (precise definitions will be given later in the text) is also derived using simple 
statistical arguments. Finally, we have a brief discussion in Section IV. 

II. THE BURST MODEL 

Consider a two-dimensional square lattice with periodic boundary conditions. (Rules for 
other lattice geometries are similar.) We denote the energy stored at each lattice site r by 
E(r)(> 0). At each timestep, the system is parallel updated according to the following rules 

1. Energy introduction: a unit of energy is added to a randomly chosen site r , i.e. 
E(t ) -> E(t ) + 1. 

2. Fire triggering: if the energy in the chosen site is greater than or equals to a fix 
threshold E c x, the energy will be dissipated (or burnt), i.e. if E(r ) > E cl , then 
E{r ) -> 0. 

3. Fire propagation: if the energy in the nearest neighbor of a burning site is greater 
than or equal to another fix threshold E c2 (< E cl ), the energy in that neighboring site 
will be also dissipated. That is, if E(r' ) > E c2 , then E(r f ) — > 0, where r' are the 
nearest neighbors of the burning site r. This process is repeated until the fire cannot 
propagate any longer. And we also assume that fire propagation is so fast that it takes 
only one timestep to set the relevant sites on fire. 
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4. Material diffusion: diffusion takes place at every site, i.e. E(r) — > E(r) + AE(r) for 
all r, where AE depends on the geometry of the lattice and the diffusion constant -DQ. 

The system is allowed to evolve for a sufficiently long time in order to attain a steady state 
before any statistics is taken. 

The behavior of the burst model depends on four external parameters, namely, E c i, E c2 , 
D, and the size of each energy blob in comparison with the size of a lattice cell. Depending 
on various choice of these parameters, the occurrence probability P(S) of dissipation size S 
falls into one of the following three types of behaviors 0: (a) P(S) shows a localized peak 
which is a signature of a super-critical state; (b) P(S) shows an early exponential cut-off 
which can be identified as a sub-critical state; and (c) power law is observed in P(S) at the 
transition point between the super-critical and sub-critical phase. 

III. STATISTICS IN THE SUPER-CRITICAL PHASE 

Let us introduce a few concepts before discussing the super-critical phase statistics and 
the relationship between the burst model at super-critical phase and percolation. A site 
is called burnable if its energy content is greater than or equals to E&. Thus, a burnable 
site will catch fire when any one of its nearest neighbors burns. Then, we define a burnable 
cluster to be a collection of maximal^ burnable sites such that all of them catch fire in the 
next timestep provided that any one of them is on fire [j/j. The size of a burnable cluster s 
is defined as the number of sites it contains. 

During the evolution of the system, clusters of sites grow continuously while they burn 
suddenly to empty spaces (or voids). And it is interesting to note that the dissipation size 
S (i.e., the total energy released in the burning of a cluster) is generally proportional to the 
cluster size. The behavior of P(S) under different choice of parameters can be understood 
in terms of the size of the burnable clusters J7j. In particular, if there is only a few large 
burnable clusters in the system immediately before most of the burnings, P(S) will show a 
peak indicating that majority of the energy release in the system occurs in the "large events" 
which extend across the whole system. Therefore, the system is in the super-critical phase. 
Finally, we define the coverage fraction p to be the percentage of sites in the system that 
found to be burnable. 

As we can see in Fig. Ill which is a snapshot of a typical system immediately after a large 
energy dissipation, a large void is observed with cluster islands of different sizes lying inside. 
So, in general, a burnable cluster for the system (in two-dimension) is not compact. 

To investigate the relation between our model and percolation, we consider a region of 
space with linear size L inside the system. If the typical length of a burnable cluster £ is 
much less than L, we expect most of the dissipation in this area is initiated and confined 
or localized inside the region. On the other hand, if £ <^ L, the coverage fraction in this 
area would gradually grow to a value pb before a large dissipation sweeps through the whole 



1 Nevertheless, our study is focused on diffusionless system in this paper. 
2 The collection of sites is maximal in the sense of inclusion. 
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region. Therefore, if the system is super-critical, we can take a sufficiently large L and 
still have non-negligible probability of finding burnable clusters with size ^ L 2 . In order to 
burn down such a large cluster, the coverage fraction around in the large burnable cluster 
immediately before it burns must be greater than or equal to the (site) percolation threshold 
p c of the corresponding lattice. In general, the (time) mean coverage fraction p depends on 
the large as well as small dissipations. But at the super-critical phase, large dissipations 
dominate the behavior of the system and hence we can safely neglect the effect of small 
dissipations. 

Using the same argument, we know that the coverage fraction p attains its local minimum 
value just after a fire. Then as energy is gradually introduced into the system, p increases. 
However, the rate of increase p will, in general, flatten out gradually with time due to the 
fact that adding energy into a site which is already burnable does not increase p. Eventually, 
p reaches its local maximum value immediately before a fire. A typical time evolution of 
p for a super-critical system is plotted in Fig. |2|. In this respect, the time evolution of 
burnable clusters in the burst model can be viewed as a dynamical percolation problem 
where the "conductance probability" gradually increases with time. Any sudden drop in the 
conductance probability is due to a fire. Finally, we like to remark that both p& and p are 
functions of E cl and E c2 . 



A. The Case When E c2 = 1 

We first consider an easy case with D = 0, E c2 = 1 and E c \ ^> E c2 . With this set of 
parameters, the system is in the super-critical phase 0. Since E c2 = 1, an empty site turns 
burnable once an energy unit is dropped onto it. Moreover, burning is the only way to turn 
a burnable site back to an empty site because D = 0. Thus, a burnable site is formed out 
of an empty site with probability (1 — p)/N at each timestep where N is the total number 
of sites in the system. That is, 

g = *(W) (i) 

for some positive constant k. Integrating Eq. (HD gives us 

Pit) = 1 - e~ kt , (2) 

where t is the time since the site is last burnt (and hence p(0) ~ 0). Since the system is in 
the super-critical phase and E c2 = 1, it is almost sure that immediately before each burning, 
almost the entire system is covered by a single large burnable cluster. So after each burning, 
almost all the sites becomes void. Everything starts over again almost independent of the 
history of the system. Thus, we may relate the mean coverage fraction, p, to the average 
coverage fraction attained just before burning, p&, by 

where the expectation is taken over all the values of coverage fractions, p', just before 
burnings. 
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In the E c i — > oo limit (while the lattice size remains fix), every grid point the system 
must be covered by a single burnable cluster right before each fire. So p b = p = 1, which 
agrees with the prediction of Eq. ([3]). As we lower the value of E c i, some sites may not be 
contained in the burnable cluster before a fire. Therefore, we expect both p& and p to be 
less than 1. At this point, Eq. (^) still holds. As we further decrease the value of E c i, the 
system eventually enters the critical, and then follows by the sub-critical phase. As there is 
negligible fire which extends system-wide, fluctuation of p about its (time) mean value p is 
small. So for a infinitely large system, p b equals p and their relation given by Eq. (§) is no 
longer valid. Finally, as E c \ — ► E c2 , both p b and p goes to 0. 

In order to study the local properties of the system near a burnable cluster, we define 
(p&) Local to be the (time) average coverage fraction immediately before a fire at a sufficiently 
small region centered at the triggering site. Using the same argument in the derivation of 
Eq. (H), we find that 

n ~ 1 Unocal ~ P(°) U) 

P ~ log(l-p(0))-log(l-( Pb ) Local ) ' 1 ) 

where p(0) is determined by the amount of energy remains after the last fire in the nearby 
region, which should be close to 0. We expect (pb) Local to be greater than its corresponding 
global value p&. While p& does not fluctuate in the limit of large system size in the sub-critical 
phase, (p&) Local near a burning cluster can deviate from p b . 

The above discussions are verified by numerical simulations^. We fix E c2 = 1 and D — 
on a 512 x 512 square lattice, p is obtained by calculating the time mean of p. On the other 
hand, p& and (pb) Local are obtained by calculating the median of their corresponding values 
immediately before each fire. Their values obtain this way will not be seriously affected by 
the "pre-mature" burnings. As shown in Table |, Eq. (|j) agrees very well with our numerical 
experiment in the super-critical phase. In addition, the relation between p and (pb) Local in 
both the super-critical and sub-critical phase is also verified. As one can notice from Table |, 
the value of p predicted from (pfe) Local deviates systematically from the observed valued of 
p in the sub-critical phase. The deviation can be explained as follows: Eq. (Q) tells us that 
the predicted value of p increases when we increase p(0) but keeping p b fixed. Consequently, 
when p is fixed, p b decreases with increasing p(0). Thus, in the sub-critical phase whose 
p(0) > 0, the prediction of p in Table | using the assumption that p(0) = is systematically 
drifted to the lower side. 

All the above observations are also verified for honeycomb, triangular, and simple cubic 
lattices (see Table |). In the supercritical phase, Eq. (|3]) holds very well in all these lattices. 

Since we only use a modest set of assumptions in deriving these relations, we expect 
Eqs. (0)-(§) to hold for a number of models with similar growth mechanisms as long as 
there is a separation of timescale between the gradual growth and sudden dissipation. To 
test this hypothesis, we look into the version of FFM proposed by Drossel et al. We 
have performed numerical experiment on a 1024 x 1024 square lattice with the probability 
of tree grow and lightning per site per timestep being 1 x 10~ 9 and 1 x 10~ 13 respectively. 



3 AU computer simulations in this paper are done in SUN SPARC 10 workstations in University 
of Hong Kong, and SGI Challenge in IAS 
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The average coverage fraction p is found to be 0.408(1) which is consistent with the value 
0.4081(7) found recently by Clar et al. 0. We also obtain p b = 0.423(2) and (p&) Local = 



0.62(3). If we assume that p(0) = 0, p deduced from the local coverage fraction equals 
0.359(2) which is again drifted to the lower side. Nevertheless, within an error of 10% or so, 
Eq. (f|) gives a reasonably good estimate of p based on (p&) Local - In this respect, we believe 
that Eq. @ is a good approximation for models with growth mechanisms similar to the 
burst model or the FFM. 



B. The Case When E c2 > 1 

The situation is more complicated when the system is in super-critical phase with E c2 > 
1. For example, if we take E c \ = 16, E c2 = 2 and D = in the 512 x 512 square lattice, 
Table |I| tells us that p b = 0.874(3) and p = 0.502(1). Although the system is in the super- 
critical phase, Eq. (^) does not hold. This discrepancy is due to the invalidity of Eq. ([!]). 
Since the adding an unit of energy to a site with E < E c2 does not necessarily turn it into 
a burnable site, we expect dp/dt to be smaller than k(l — p). 

For simplicity, we consider a diffusionless system and write E c2 — n - We also define pj to 
be the coverage fraction of the system by sites with energy equals j for j = 0, 1, 2, . . . , n — 1. 
Moreover, p n is the coverage fraction of the system by sites with energy greater than or 
equals to n (hence they are burnable and p n = p). It is easy to see that 

IF = ~ kp0 ' (5a) 

^fl = kpj_! - kpj for j = 1, 2, . . . , n - 1, (5b) 

and 

dp n 
dt 

where k is a positive constant. Note that the coverage fractions satisfy the conversation law 



kpn-i , (5c) 



EP, = 1- (6) 

Clearly, Eq. (||c) reduces to Eq. (]l|) when n = 2. 

Given the initial coverage fractions, we can integrate Eqs. (|S|) one at a time, which can 
eventually calculate the time evolution of pj for all j < n. Then p n (t) can be computed 
using Eq. (J^j. After some simple but long computation, the solutions of Eqs. (0) are given 
by 

3 (kt)' m 

Pi® = E M-Pi-m (0)e- fe * for j = 0, 1, 2, . . . , n - 1, (7a) 



m=0 m - 



and 
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n-1 3 (L,f\m 

pn(t) = i - E E ^rP,- m (o)e" fct , (7b) 

where t is the time since the last burning. Just after a large energy dissipation, po(0) ~ Pb, 
Pi(0) ~ (1 — Pb)/( n — 1) for j = 1, 2, . . . , n — 1, and p n (0) = are reasonable approximations. 
Thus, Eq. (^|b) becomes 

Pn(t) « 1 - * £ -^£^(»-2- j)e~ kt • (8) 

3=0 J" 3=1 

In particular, if E7 cl is very large, p fe w 1, and Eq. becomes 

p„(t) - 1 - e-- £ ^ . (9) 

m=0 

Unfortunately, there are no simple close forms for the time T when p(T) = p n (T) = 
in general for both Eqs. @ and (Q). So we do not have a simple relationship between p and 
as in the case when E c2 = n = 1. 

To verify the validity of Eq. @, we numerically solve for fcT when p n (T) = p^, and 
n = 2. Once fcT is found, we can numerically integrate Eq. (|8p to obtain p. The results are 
tabulated in Table ||. As we can see, the prediction of Eq. (H) agrees quite well with our 
numerical experiments. 



IV. DISCUSSIONS 

In summary, we have investigated the statistics of the burst model at its super-critical 
phase. We have a simple picture for the time evolution of the coverage fraction p. Using 
simple statistical argument, we also derive a relationship between the (time) mean coverage 
fraction p and the average coverage fraction just before a fire pb- The relation agrees quite 
well with our numerical simulations within an error level of a few percent. Moreover, we find 
a local version of this relationship (see Eqs. ©-©)• Since these two relations are derived 
using a modest set of assumptions, we expect that they are valid in many other growth 
models with dissipations. Our preliminary analysis on the FFM tells us that this is indeed 
the case. And it will be very interesting to apply the same analysis to both the burst model 
in higher spatial dimensions and other growth and dissipation models. 
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TABLES 



TABLE I. Values 


of the averaj 


$e coverage 


tractions 


when D = (J 


and h C 2 = 


1 

1 m various 


underlying 


lattices with different values of E c \ . 










Lattice 


Ed 


/) 

r 


Ol, 


VPWLocal 


/ — \ a 
Wprojl 


/ — \ b 


Phase 


512 X 512 


4 


0.334(2) 


0.338(4) 


0.55(4) 




0.31(2) 




Square 


6 


0.381(1) 




0.63(3) 




0.36(2) 


Sub-critical 




8 


0.400(4) 


0.41(1) 


0.66(1) 




0.388(7) 


T 




10 


0.426(5) 


0.61(1) 


0.71(1) 


0.35(1) 


0.426(8) 


I 




11 


0.509(5) 


0.79(1) 


0.800(4) 


0.49(1) 


0.503(4) 


Super-critical 




12 


0.580(4) 


0.872(5) 


0.872(4) 


0.576(6) 


0.576(5) 






13 


0.633(3) 


0.920(3) 


0.920(3) 


0.636(4) 


0.636(4) 




cin v c-i o 
01Z X 0±Z 


8 


0.473(2) 


0.476(2) 


0.73(4) 




0.44(3) 


Sub-critical 


Honeycomb 


10 


0.495(4) 


0.52(1) 


0.76(3) 




0.47(3) 


T 




12 


0.530(3) 


0.81(2) 


0.82(1) 


0.51(2) 


0.52(1) 


I 




14 


0.666(1) 


0.942(2) 


0.941(1) 


0.669(3) 


0.668(2) 


super-critical 


512 X 512 


6 


0.323(4) 


0.326(4) 


0.52(4) 




0.29(3) 


sub-critical 


Triangular 


8 


0.337(3) 


0.400(5) 


0.54(3) 




0.30(3) 


T 




10 


0.465(1) 


0.753(4) 


0.754(4) 


0.462(3) 


0.462(3) 


I 




12 


0.586(1) 


0.880(2) 


0.880(3) 


0.584(2) 


0.584(3) 


super-critical 


64 x 64 x 64 


6 


0.220(2) 


0.230(6) 


0.30(3) 




0.15(2) 


sub-critical 


Simple Cubic 


8 


0.294(4) 


0.488(4) 


0.510(5) 


0.271(3) 


0.285(4) 


I 




10 


0.464(2) 


0.756(2) 


0.756(4) 


0.464(2) 


0.464(3) 


super-critical 



a (<°)proji ^ S t ne value of p obtained as if Eq. (||) is correct. The difference between this value and 
p we measure (in column 3) tells us how good Eq. (|3|) is. We only show this prediction when the 



system is in super-critical state. 

b (p) pro j 2 is the value of p obtained by substituting (/0fe) Loca i into Eq. (Q) and assuming p(0) = 0. 
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TABLE II. Values of the average coverage fractions when D = and E& = 2 in 512 x 512 
square lattice using different values of E c \ . 





P 


Pb 


\Pb) Local 


(p)pro]l 8 ' 


(/ 3 )proj2 b 


Phase 


8 


0.338(1) 


0.338(1) 


0.60(5) 




0.32(4) 


sub-critical 


10 


0.358(2) 


0.361(3) 


0.64(4) 




0.34(4) 


T 


12 


0.368(2) 


0.387(3) 


0.68(1) 




0.365(8) 


I 


14 


0.379(2) 


0.557(3) 


0.685(6) 


0.297(3) 


0.368(5) 


super-critical 


16 


0.502(1) 


0.874(3) 


0.875(5) 


0.503(3) 


0.504(5) 





a ( / o) pro j 1 is the value of p obtained as if Eq. (||) is correct. The difference between this value and 
p we measure (in column 2) tells us how good Eq. (|8|) is. We only show this prediction when the 
system is in super-critical state. 

b (p) P roj2 i s * ne va l ue °f P obtained using the local version of Eq. @. 
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FIGURES 



FIG. 1. Gray scale snapshot of a 512 x 512 system with E c \ = 8, E C 2 = 1 and D = just after 
a large dissipation showing a large void inside which there are some small burnable clusters. 

FIG. 2. The time evolution of the coverage fraction p (in solid line) in a typical super-critical 
system. Here, we use E c \ = 13, E C 2 = 1, and D = 0ina512x512 square lattice. The values of p 
and pb are shown (in dash lines) for comparison. 
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